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1. Introduction 



One of the long standing questions in astrophysics is the progenitors of Gamma Ray 
Bursts (GRBs). Collapse of a massive star is one of the most favoured progenitors of the 
long GRBs. Evidence for a massive star being a GRB progenitor may be obtained in two 
different ways, both using observations of GRB afterglows. (1) Supernova ( SN) component 
underlying the GRB afterglow : A few of the nearby GRB afterglo ws have shown the tem- 
poral and spectroscopic signature of an underlying SN. see e.g., (jStanek et al.l 120031 ) (2) 
Evolution of GRB afterglow in the stellar wind medium : Massive stars modify the density 
profile of the circum-stellar medium due to the powerful winds they drive during their life 
time. For a constant mass loss rate and constant wind velocity the circum-burst medium 
assumes a density profile p oc r~ 2 as compared to p = constant in the absence of stellar wind. 
Evolution of th e GRB afterglow light curves is significantly diffe rent in these two cases of 
density profiles (jWijers and Galamalll999l ; IChevalier and Lill2000l ). Attempts to look for the 
signatures of such a wind-modified circum-burst density profile in the light curves o f GRB 
afterglows have not been conclusive so far. In the case of GRB 050904 iGendre et al.l (120071 ) 
finds that the early x-ray afterglow suggests a wind like density profile of the circum-burst 
medium while the late optical afterglow was consistent with evolution in a constant density 
medium. Hence, they conjecture that a transition between these two types of density profiles 
would have taken place somewhere in between. However, this transition was not directly 
observed in the light curve of any given band. We show that the afterglow of GRB 050319 
could be explained as being due to the transition of the circum-burst density profile from 
wind-like to constant density. 



GRB 050319 was detected by the Burst Alert T elescope (BAT) inst rument of the 
SWIFT satellite on 2005 March 19, 09:31:18.44 UT fcrimm et al.l l2005al lbh. However, 
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Cusumano et al.l (120061 ). using the re-analysis of the BAT data, pointed out that Swift 



was slewing during the GRB onset and the BAT trigger was switched off. The GRB was 
recognised about 135 seconds after its a ctual onset. The tota l duration of the GR B (T 90 ) 
was thus 149.7 s (jCusumano et al.l 120061 ) instead of 10 ± 2 s (jKrimm et al.ll2005al lbT). The 
burst fluence in 15-350 keV band within the T 90 duration is estimated to be 1.6 x 10~ 6 
erg cm -2 . The photon index of the time-averaged single power law spectrum is 2.1 ± 
0.2. Swift XRT and UVOT located a bright source at a = 10 h 16 m 48 s and 8 = +43°32 W 



(J2000) which was later confirmed by iRvkoff et al.l (120051 ) with ROTSE-IIIb. iFvnbo et al. 



(120051 ) obtained the spectra of the afterglow of GRB 050319 on 2005 March 20 and the 
redshift of the afterglow was measured to be z = 3.24. At this redshift, the gamma ray 
isotropic equivalent energy released during the burst is 3.7 x 10 52 erg for a flat universe with 
Vl m = 0.3, Vl A = 0.7 and H = 70 km s" 1 Mpc" 1 . 



Observations and Data Reduction 



Optical CCD observations of the afterglow of GRB 050319 were carried out in Johnson 
BV and Cousins RI filters using the 104-cm Sampurnanand Telescope of ARIES, Nainital 
with regular specifications of the CCD camera an d using standard ob servation procedures 
of bias subtraction and flat fielding. For details see iMisra et al.l (120071 ). 



The BVRI magnitudes of the optical transient (OT) obtained were c alibrated differ- 
entially using secondary stars numbered 8, 9, 10, 11 and 13 in the list of iHendenl (120051 ). 
The magnitudes derived in this way are given i n Table ffl Photome t ric magnitudes avai l- 
able in the l i teratu re by IWozniak et al.l (120051 ) ; iQuimby et ajj (120061 ) ; iMason et al.l (120061 ) ; 
Huang et al.l (120071 ) were converted to the present photometric scales using the five secondary 
stars mentioned above. 



3. Light curves of GRB 050319 



Along with our own observations we have used observations reported in the literature to 
study the light curves of GRB 050319 af terglow. The x-ray afte rglow was observed by Swift 
XRT starting from ~ 220s to 28 days (jCusumano et al.l 120061 ) after t he burst. Attemp ts 



to observe the afterglow at r adio wavebands resulted in upper limi ts dSoderberg 



2005a 



Volvach and Pozanenko 2005). The optica l afterglow was observed by Wozniak et al. " j2005 ): 
Quimby et al.l (120061 ); IMason et al.l (120061 ) resulting in a coverage from a few seconds to ~ 4 
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Table 1: CCD BVRI broad band optical observations of the GRB 050319 afterglow using the 
104-cm Sampurnanand Telescope at ARIES, Nainital. 



days after the burst. 

To construct the optical light curve w e have correct ed the observed magnitudes for the 
standard Galactic extinction law given by iMathis! (119901 ). The galactic extinction in the di- 
rection of GRB 50319 is estimated to b e E(B-V) = 0.011 mag from the smoothed reddening 
m ap provided by ISchleeel et al.l (119981 ). The effective wavelength and normalization given 
by Bessell et al. ( 19981 ) were used to convert the magnitudes to fluxes in /iJy 



Most of the GRB afterglow light curves are well characterized by a broken power law 
of the form F = F {(t/tb) QlS + (t/t b ) a2S } _1 / s where a,\ and a 2 are the afterglow flux decay 
indices before and after the break time (tb), respectively. Fo is the flux normalization and 
's' is a smoothening parameter which controls the sharpness of the break. Most known GRB 
afterglows have a 2 > «i i.e. the decay becomes steeper after the break. Interestingly, the 
optical afterglow light curve of GRB 050319 shows steeper to flatter decay with a break at 
~ 0.02 day. This behavior of light curve decay is difficult to explain within the standard af- 
terglow models. The x-ray and optical light curves also show some variability superimposed 
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on the power law decay. The x-ray light curve shows a break near ~ 0.3 day. We quantify 
the various characteristics of the afterglow light curves as summarized below. 



1. The x-ray afterglow of GRB 050319 shows a very rapid decay before ~ 0.005 day 
^ 384 s) after the GRB . The decay then flattens before steepening again at ~ 0.3 day. 



Cusumano et al.l (120061 ) has characterised the afterglow into three separate temporal 
evolutions in x-ray bands : a X i = 5.53 ±0.67 (At < 384 s); a X2 = 0.54 ±0.04 (384 s < 
At < 0.3day); a X3 = 1.14 ± 0.2 (At > 0.3day); 

2. The B, V and R band light curve also show a rapid decline during the early phase 
(At < 0.02 day) which flattens at later epochs. Thus, the afterglow can be separated 
into two separate temporal evolutions in R band ori = 1.09 ± 0.03 (At < 0.02day); 
a R2 = 0.51 ± 0.03 (0.02day < At < lO.Oday); 

And for B, V and I bands we measure following a s : q;bi = 1.46±0.26 (At < 0.02day); 
a B2 = 0.33 ± 0.05 (0.02day < At < l.Oday); a V i = 0.90 ± 0.05 (At < 0.02day); 
a Y2 = 0.50±0.04 (0.02day < At < l.Oday); a l2 = 0.59±0.13 (0.02day < At < l.Oday); 
The average decay index of the optical light curve, at early and late epochs, is then 
1.15 ± 0.27 and 0.48 ± 0.15 respectively. 



4. GRB 050319 afterglow : wind or homogeneous density profile ? 

The breaks seen in x-ray light curve (at ~ 384 s and at ~ 0.3 day) are not accompanied 
by simultaneous breaks in optical wavebands. Similarly, the break seen in optical band has 
no simultaneous counterpart in the x-ray light curve. Also, the sense of slope change, i.e. 
a 2 < Qi, as see n in o ptical waveband is contrary to the predictions of the fireball model 
rtSari et aliE998l . ll996h which expects at\ < a 2 . 



We propose a different model to explain the afterglow of GRB 050319 as being due to 
a transition of the circum-burst medium density profile that the explosion generated shock 
wave is interacting with. We propose that the observed change from steep to flat decay of 
the optical afterglow of GRB 050319 at 0.02 day is due to the change in the density profile 
of the circum-burst medium from wind modified (p oc r~ 2 ) to the constant density medium 
(p = constant). The break in the light curve occurs when the shock front interacts with 
the boundary between the two density profiles. Below we describe and reproduce various 
features of the GRB 050319 afterglow using this model of 'wind to constant density medium 
transition'. 
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The early steep decay of x-ray afterglows (a ~ 3 to 5) as is seen i n the case of GRB 
050319 are now seen routinely in most of the GRBs ( INousek et al.ll2006l ) a nd has become a 



canon ical feature of the GRB x-ray afterglows. In the case of GRB 050319 ICusumano et al. 



( 120061 ) conjecture that the early steep decay emission could be low energy tail of the GRB 
prompt emission. We exclude this early emission from the rest of our discussion and we will 
restrict ourselves to the rest of the x-ray light curve. 

The radiation spectrum of GRB afterglows exhibits a power law spectrum characterised 
by three break frequencies - the self absorption frequency z/ a , the peak frequency v m corre- 
sponding to the lower cutoff in the electron energy distribution (71(7) oc 7~ p , 7 > j m ), and 
the synchr otron cooling fr equency u c . The flux F m at u m provides the normalisation of the 
spectrum ( ISari et al.lll998l ). 



The photon index (r) of the afterglow and the electron energy distribution index p 
are related in any given spectral regime (r — 1 = p/2 if v c < v and T — 1 = (p — l)/2 if 
v < u c ). The corresponding temporal decay index a would be (3p — 2)/4 and 3(p — l)/4 
respectively before the jet break and would equal p in both spectral regimes after the jet 
break, according to the standard fireball model for an afterglow expanding in a homogeneous 
interstellar medium (ISM). For the shock wave expanding into the wind density profile the 
corresponding a would be (3p — 2)/4 and (3p — l)/4 respectively before the jet break and p 
after the jet break. 

In the present case, the observed values of the photon index (T = 1.69 ± 0.06) and 
temporal decay index (ax = 0.54 ± 0.04) of the x-ray afterglow are consistent with the 
spectral regime ux > v c and p = 1.5. The observed decay indices of optical light curve 
are also consistent with the inferred value of p, the spectral regime i/ pt < and wind- 
constant density transition at 0.02 day. As discussed above, the expected temporal decay 
index of the x-ray afterglow, «x — (3p — 2)/4, is the same for wind and homogeneous ISM 
density profiles. The absence of a break in the x-ray afterglow light curve simultaneous 
with the optical break makes the multiband afterglow features consistent with the proposed 
transition of the circum-burst medium density profile from wind to constant density. In 
Figure 1. we compare the predictions of our model of wind to constant density transition 
of the circum-burst density profile with multiband observations of GRB 050319 afterglow. 
A detailed list of the best fit spectral parameters can be found in Table [2j The observed 
B band light curve is systematically lower than that predicted by the model which could be 
du e to the Lym a n a a bsorption at z = 3.24 appearing in the observer's B band as suggested 
bv lHuang et al.l (120071). The steepen ing of the x-ray light curve at At ~ 0.3 day could be due 
to jet break (ICusumano et al.l 120061 ) but unfortunately the variability in R band light curve 
and insufficient sampling of the data in B and V bands after ~ 1.0 day makes it difficult to 
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verify the achromaticity of the break. 

Given the above model spectral parameters we find that the broadband behavior of the 
afterglow is very well explained. However, we are restricting ourselves to the overall behavior 
of the afterglow and hence do not attempt in our model to reproduce the variations seen in 
the optical light curve. The reason for these variations could be density inhomogenities in 
the circum-burst medium. 



4.1. Physical Parameters 



Four spectral parameters (z/ a , u m , u c and F pea k) are related to four physical parameters 
viz n (no density of the constant density circum-burst medium) or A* (defined as p(r) = 
5 x 10 11 A^r -2 for wind density medium), E (total energy content of the fireball), energy 
fraction in relativistic electrons e e and that in magnetic field e#- The typical value of self 
absorption frequency u a lies in radio-mm waves and hence is best estimated only if the 
afterglow is well observed in these ba nds. Unfortunate ly, the afterglow of GRB 05031 9 
was never detected at the radio band (jSoderbergJ l2005al ibl; IVolvach and Pozanenkol 120051 ). 
Therefore we expressed the remaining three spectral parameters, known separately from parts 
of the light curves corresponding to the wind and the constant density circum-burst medium, 
in terms of A* and n, respectively. Equating the kinetic energies estimated from two density 



profiles, i.e. E, 



K 

wind 



& 'ISM 'i 



we obtain a relation between A* and n : A* = 5.097 x 10" 



n 



2/5 



For a typical range of values of n (1 to 100), estimated A* ranges from 5.097 x 10~ 3 to 0.032 
which results in the range of Ef so from 1.3 x 10 54 to 5.3 x 10 53 erg. All the estimated physical 
parameters are listed in Table [3J 



5. Discussion 



5.1. Signature of Wind Reverse Shock ? 



Morphology an d evolution of the wind bubbles has been studied by lCastor et al.l (119751 ); 
Weaver et al.l (119771 ). A reverse shock forms at the surface where the stellar wind meets the 
surrounding ISM and it then propagates into the wind. The free wind (upstream of the 
reverse shock) has a density profile p oc r -2 and the shocked wind (downstream of the reverse 
shock) has a constant density pr ofile. The effects of such a density transition on the afterglow 
of a GRB has been studied by iPe'er and Wijersl (120061 ). It could be this transition of the 
density profile that we are observing at 0.02 day in the present case of GRB 050319. From 
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the observations of long GRB afterglows, it has been inferred that most of the GRBs occur 
in constant density environment and the absence of wind signatures in the GRB afterglow 
was surprising. Various ways which can bring the wind reverse shock closer to the exploding 
star have recently been proposed to r esolve this mystery surrounding the absence of winds 
( Ivan Marie et al.ll2006l ; lEldridgd 120071 ). In the case of GRB 050319, for a range of assumed 
values of n = 1 to 100 we estimate the radi us of the r everse shock to be Rsw ~ 0.5 pc to 
0.1 pc comparable to the values obtained by lEldridgd (120061 ). The constraint t B < 1 puts a 
lower bound on density : n > 0.03 



5.2. Implications for the models of GRB progenitors 



Our interpretation of the afterglow of GRB 050319 as being due to the wind-constant 
density transition supports the collapsar model of GRBs. Detection of similar transitions in 
other afterglows have so far eluded us perhaps because of the smaller size of the wind bubbles 
and the resultant early transition times. The present detection was ma de possible chiefly be 
cause of the quick follow up abili 



and ROTSE-IH (IQuimbv et al 



; ies of the robotic telescopes RAPTOR (IWozniak et al. 



2006|) and UVOT (IMason et al 



2005) 



2006) coupled with those of Swift XRT (jCusumano et al. 
2006h . Time dilation due to cosmological redshift delays the 



occurrence of the transition in the observer's frame of reference and makes it favourable to 
detect such a transition in distant GRBs. The robotic telescopes are now routinely detect- 
ing GRB afterglows as early as a few minutes after the burst and with careful analysis of 
multiband observations of distant GRBs it should be possible to detect more examples of 
similar transition. For example, a probable detec tion of similar densi ty transition, though 
not in the same waveband, has been reported by iGendre et al.l (120071 ) in the case of GRB 
050904 (z ~ 6.3) where the transition time is assumed to be ~ 1700 s (~ 0.02 days) after 
the burst, similar to that for GRB 050319 afterglow in the present case. 



6. Summary 

We have modeled the multiband afterglow of GRB 050319, using our own optical ob- 
servations and other observations available in the literature, as being due to the interaction 
of the relativistic blast wave with circum-burst medium which shows a transition of density 
profile from wind to constant density. Our conclusions can be summarised as follows : 

1. We present BVRI band observations of GRB 050319 afterglow. 

2. We showed that the unusual break in the light curves of optical afterglow at 0.02 
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day can be explained as being due to the transition of circum-burst density profile 
from wind to constant density. The observed x-ray afterglow light curve without a 
simultaneous break is consistent with this interpretation. The overall afterglow can be 
explained by using a relatively low value of electron energy distribution index p which 
is also consistent with the x-ray spectral photon index. 

3. The transition of the density profile could be due to the wind reverse shock propagating 
into the stellar wind driven by the progenitor of GRB 050319. We estimate radius of 
the wind reverse shock to be Rsw ~ 0.5 pc to 0.1 pc for assumed values of n ~ 1 to 100 
cm -3 respectively. 



Acknowledgments 

We are thankful to D. Bhattacharya for critical comments and detailed discussions 
throughout this work and to Ram Sagar for his support during the observations. We are 
thankful to the anonymous referee for constructive comments which has improved the Letter 
significantly. 

REFERENCES 

K. Z. Stanek, T. Matheson, P. M. Garnavich, P. Martini, P. Berlind, N. Caldwell, P. Chal- 
lis, W. R. Brown, R. Schild, K. Krisciunas, et al, ApJ Lett 591, L17 (2003), 
astro- ph/0304173. 

R. A. M. J. Wijers and T. J. Galama, ApJ 523, 177 (1999), |astro-ph/9805341 . 





wind density medium 


constant density medium 


v m {Hz) 


1.61°;? x 10 13 


1.0+°; 35 x 10 12 




l.l+ll x 10 15 


2.1±g-| x 10 15 


Fpeaki^Jy) 


2370 ± 355 


18101??° 


p 


1.59 ±0.06 


1.52 ±0.02 


X 2 dof( d °f) 


1.4 


t (161) 



Table 2: The best fit spectral parameters for wind (< 0.02 day) and constant density (> 0.02 
day) profile. All the parameters are fitted at 0.003 day after the burst. 
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constant density medium at 0.02 day. The best fit spectral parameters of this model are 
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